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ABSTRACT 

In four globular clusters (GCs) a non negligible fraction of stars can be in- 
terpreted only as a very helium rich population. The evidence comes from the 
presence of a "blue" main sequence in oj Cen and NGC 2808, and from the the 
very peculiar horizontal branch morphology in NGC 6441 and NGC 6388. Al- 
though a general consensus is emerging on the fact that self-enrichment is a com- 
mon feature among GCs, the helium content required for these stars is Y ^ 0.35, 
and it is difficult to understand how it can be produced without any — or, for 
oj Cen, without a considerable — associated metal enhancement. We examine the 
possible role of super-AGB stars, and show that they may provide the required 
high helium. However, the ejecta of the most massive super-AGBs show a global 
CNO enrichment by a factor of ~4, due to the dredge-out process occurring at 
the second dredge up stage. If these clusters show no evidence for this CNO 
enrichment, we can rule out that at least the most massive super-AGBs evolve 
into O-Ne white dwarfs and take part in the formation of the second generation 
stars. This latter hypothesis may help to explain the high number of neutron 
stars present in GCs. The most massive super-AGBs would in fact evolve into 
electron-capture supernovae. Their envelopes would be easily ejected out of the 
cluster, but the remnant neutron stars remain into the clusters, thanks to their 
small supernova natal kicks. 

Subject headings: globular clusters: general — stars: formation — stars: AGB 
and post-AGB — stars: neutron 



-2- 



Introduction 



The majority of t he inhom o geneit i es in the chemi c al com position of Globular Cluster 



(GC) stars — see e.g. ISneden I (119991 ); iGratton et al.l (120041 ) — appear due to primordial 



enrichment by hot-CNO cycled material processed in stars belonging to a fi rst stellar gen- 



eration. Either massive AGB envelopes subj ect to hot bottom burning, (e.g. IVentura et al. 



(FRMS) ( 



200ll ; I Ventura. D'Antona, fe Mazzitelrill2002r) or the envelopes of massive fastly rotating stars 



Decressin et al. 20071 ) are an ideal environment to manufacture elements through 



nuclear reactions in which proton captures are involved. A second stellar generation would 
then be born from the ejecta of these stars, either directly, or diluted with pristine material. 
There are however problems with these scenari o. For the AGB progenitors, problems con- 
cern t he chemistry of the anomalous stars (see Denissenkov fc HerwieJ 12003 ; Ventura et al. 



2004 ; iFenner et al.l 120041 . but consider also I Ventura &; D'Antonal 120051 ) , and the require- 
ments for the mass budget necessary to pro duce a second stellar generation that, today, is 



generally about as numerous as the first one (ID'Antona &; Caloill2004l ; iBekki &; Norrisl 12006 



Prantzos fc Charbonnelll2006l ). The FRMS model is somewhat less problematic concerning 
the constraints on the IMF, but poses, in any case, severe problems for "normal" GCs, that 
show a unique metallicity for all stars (with the notable exception of u Cen): this model 
requires that the envelopes of these stars contribute to the second star formation stage, but 
the supernova ejecta are expelled from the clusters and do not enrich the second generation 
matter by heavy elements. 

Both scenarios predict that the stars showing chemical anomalies also must have en- 
hanced helium abundance. Helium and the hot-CNO products in the FRMS come to the stel- 
lar surface by means of the chemical mixing associated with the transport of angular momen- 
tum through the stellar radiative layers, during the phase of hydrogen burning. In the case of 
AGBs, the helium enrichment is mainly due to the second dredge up (2DU) which precedes 
the thermal pulse phase. Also the third dredge up (3DU) contributes to helium enrichment 
in the envelopes, but a further observational constraint is that th e sum of CNO abundance s 
remains remarkably constant in normal and anomalous stars (e.g. lCohen fc Melendezll2005l ). 
Consequently, the masses which can be responsible for the chemical enrichment must be high 
enough to be subject to a few episodes of the third dredge up, that brings into the envelope 
the products of helium burning. 

A variation in helium content is immediately reflected on the morphology of the horizon- 
tal branch (HB), which amplifies any evolutionary difference among the cluster stars. Helium 
enrichment by a small factor (up to Y~ 0.28 — 0.3 0) allowed in fact an easy interpretation of 
some puzzles posed by the HB (blue tails, gaps) (ID'Antona et al.l 120021 ; ID'Antona fc Caloi 
20041 ) . Helium variations from the MS location are much harder to be detected, but the 
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clustering of ~20% of u Cen stars into a "blue" MS (IBedin et al.ll2004r) was soon interpreted 



as a MS having an abnormally huge helium content (Y~0.40) (INorrisI 120041; iPiotto et al. 
2005h . In contrast to u Cen, NGC 2808 does not show a metallicity spread ( Carretta et al. 



20061). Nevertheless , evid ence for the presence of a blue MS in this cluster has been found 
by iD'Antona et al.l (120051 ). who interpreted it, again, as built up by stars having Y~0. 40, 
and were able to correlate them with the bluest side of the HB. IPiotto et al.l (120071 ) showed 
that the MS of this cluster is actually made up by three separate MSs, the bluest of which, 
including 10-15% of stars, must have a helium content Y~ 0.35 — 0.40. The presence of blue 
MS has not yet been found in other clusters: observationally it is a very difficult problem, 
and, in addition, this may be a characteristic of the most massive clusters only. 

Huge helium overabundances are also needed to explain all the m ain features of the 
anom alous HB of the metal rich clusters NGC 6441 and NGC 6388 JCaloi fc D'Antona 
20071 ) including: 1) the long extension in luminosity of the red clump; 2) the fact that RR 
Lyrs have a very long average period, which is unusual for a cluster of high metallicity; 3) 
the extension into the blue of the HB. This analysis also requires that ~ 15% of the HB 
population has Y>0.35. 

All these observations and their interpretation in terms of helium enrichment must be 
taken very seriously, but no chemical enrichment mechanism is able to produce the huge 
amount of helium required for about 15%-20% of the stars in these four clusters with- 



out d ramatically impacting upon the metal abundance (IKarakas et al.ll2006t iRomano et al. 



20071 ). It is clear that the presence of this extreme population requires some very particu- 
lar dynamical and chemical conditions to be understood. Nevertheless, the first step is to 
look for progenitors which may at least provide the required high helium abundance in their 
ejecta. 



The most massive FRMS may b e able to eject huge qu antities of helium ( iDecressin et al. 



20071 ). Concerning the AGB stars, IKarakas et al.l (120061 ) pointed out that Y~0.40 is not 
predicted by AGB models, and this could, in the end, rule out AGB stars as progenitors 
of the second stellar generation. A sc arne Y^0.29 w a s the maximum abundance computed 
from AGB evolution i n the models by IVentura et al.l (120021 ) for 5.5M^, while a s i ngle y alue 
Y=0.375 is quoted bv IKarakas et al.l J2OO6J) for a 5M model bv ICampbell et al.l (l2004h . P. 
Ventura & F. D'Antona (in preparation) have now completed new computation of massive 
AGBs, whose results revise upward the previous ones (see later). The maximum AGB 
mass in these recent computation is 6.3M , and the possible role of super-AGB stars has 
not yet been considered. Such models are difficult to be computed, as Carbon burning in 
semi-degeneracy has to be followed, and the uncertainty in the mass loss rate hampers any 
conclusions on their final evolution, either towards electron-capture supernova (ecSN) or 
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towards O-Ne white dwarfs. 



In the light of recent super-AGB stars models described in iPumol (120061 ). in this paper 
we discuss the possible role of these stars for the self-enrichment of GCs. We discuss both 
the helium abundance and the total CNO abundance versus core mass which can be built up 
by the most recent AGB and super-AGB models. We find out that the matter lost by the 
envelopes of super-AGBs may indeed be a good candidate for the high helium population 
of the four GCs discussed above. However, in the higher mass tail of the stars developing a 
O-Ne core the total CNO is increased by a factor ~ 4. Observations of the super He-rich 
stars may be able to check this result. 

If all the stars in the GCs having a high-He population have standard CNO, as found for 
all the chemically anomalous stars in many GCs, the most massive super-AGBs must evolve 
into ecSNe. This conclusion would be of further support to the idea that the presence of large 
population of neutron stars in GCs implies that they were born with super nova-kicks low 



enoug h not to be ejected by the cluster, as expected by the ecSN event (e.g. Ilvanova et al. 



2003). 



2. Super-AGB stellar evolution 



AGB stars develop a degenerate carbon-oxygen core and mass loss prevents them to 
reach the mass limit for degenerate explosive carbon ignition. Super-AGBs are defined 
as stars which ignite carbon in conditions of semi-degeneracy, thus non explosively, but 
are not able to ignite hydrostatic Neon burning in the resulting O-Ne core. Co nsequently, 



degeneracy increase s in the core and t hese stars may undergo thermal pulses (e.g. Ilben et al. 



19971 ; iRitossa et al.l Il999i ; ISiesa 120061 ) and lose mass as "normal" (but quite massive and 
luminous) AGB stars. Depending on the competition between the m ass loss rate a nd the 
core growth, they may then evolve into massive O-Ne white dwarfs ( INomotolll984l ) or as 
ecSN, electrons being captured first by 24 Mg and 24 Na then by 20 Ne and 20 F nuclei when the 
core mass reaches the Chandrasekhar mass. The 2DU occurs also in s uper-AGBs and, as i n 



Pumo fc Siess 



2006), 



the massive AGBs, it reduces the mass of the H-exhausted core (e.g. 
while the helium abundance increases in the external hydrogen envelope (e.g. ISiess fc Pumo 
20061 ). If mass loss is strong enough to avoid the super-AGB evolving to ecSN, as soon as 
super-AGBs evolve in the GCs, the supernova epoch is finished, and we may hypothesize 
that the second phase of star formation begins from the super-AGB ejecta. 
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3. Helium and CNO in massive AGBs and super— AGBs 



As a function of the Carbon Oxygen core mass M c , we report in Figure [T] the helium 
content in the AGB ejecta of metallicity Z=10~ 3 , for the models from 4 to 6.3M by P. 
Ventura & F. D'Antona (in preparation) including q-el ements enhancement i n the opacities, 
and core overshooting. An important difference from IVentura et al.l (120021 ) models is the 
following: in the present models, also overshooting at the bottom of the convective envelope 
is included, adopting the same parameter which describes core overshooting. This latter 
is calibrated to account for the width of the MS of intermediate mass stars. The effect of 
the envelope overshooting is to enhance the extent of the 2DU: fixed the initial mass, the 
remnant core mass at the beginning of the AGB evolution is then smaller, and the resulting 
helium content in the envelope is lar ger, as can be seen by comparing the values of Figure [1] 
with Table 1 in IVentura et al.l (120021 ). 



We also plot the helium co ntent at the 2DU for the super-AGB models of metallicity 
Z=10 -3 from 7.5M to 9.5M (jPumol 120061 ) computed without overshooting. In this case 
M c is the C-0 mass after the 2DU. M c is indeed the important physical parameter which 
determines the final fate of the star. The initial mass difference between the AGB models with 
overshooting and super-AGB models without overshooting is, for the same core mass, ~1M . 
For super-AGBs, we do not have the helium mass fraction of the ejecta, but these stars do 
not go through efficient episodes of 3DU, so the plotted value is a good approximation. As 
we see, the super-AGB models reach a helium content in the ejecta larger than the standard 
AGBs, and well approaching the values needed to be consistent with the super-He rich stars 
in GCs. The new AGB models an d the super-AGB models show n here, computed by a 
different program and input physics (jSiessll2006l ; ISiess &: Pumoll2006l ) are in good agreement 
with each other, at the same core mass. 

From Figure [1] we see that indeed the super-AGBs could be the progenitors of the 
very high helium population found in the most massive GCs. For a Salpeter's initial mass 
function, the mass budget in the ejecta from 6.5M (the minimum mass for super-AGB 
evolution, considering models with overshooting) to 8M (see later) is about 50% of the 
mass budget contained in the ejecta of normal AGBs from 5M to 6.5M . Therefore, from 
the super-AGBs mass budget, a 15% fraction of very high helium stars can be born, with 
the same mechanism which c an give origin to anoth er ~30% of moderately helium rich stars, 
as in the cluster NGC 2808 (ID'Antona et al.ll2005l ). How these stars can form directly from 
the super-AGB ejecta, so that their helium abundance remains as large as the required 
Y ^ 0.35, is a different problem. 

Figure CD also shows the total CNO enrichment F C no, the sum of CNO abundances with 
respect to their initial sum, in the ejecta of AGBs and at the 2DU for the super-AGBs. The 
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AGB enrichment is larger than ~2 for M C <O.9M , that is for M<5M . As we have noticed, 
the sum of CNO abundances apparently is constant (within a factor 2) in GC stars. On this 
basis, we should conclude that stars of M<5M should not contribute to the second stellar 
generation -unless their ejecta are diluted with pristine matter, so that the constancy of CNO 
is preserved. For the super-AGBs, we also find a sharp increase of Fqno, reaching a factor 
4 for the largest mass in exam. This is due to the occurr ence of the so called dredge-out 
process in this model (jlben et al.lll997l ; iRitossa et al.lll999l ). The helium mass fraction and 
the overall CNO content are the most robust predictions from these theoretical computations, 
due to the large uncertainties associated to both the extension o f the 3DU and the strength 



of hot bottom burn ing; (HBB) within super-AGB models (e.g. Ilzzard fc Poelarendd 12006 
Pumo fc Siesd 120061 ) . We also comment the individual abundances providing Fcno- The 
CNO abundances given for M < 6.3M are the average mass fractions of the ejecta. As HBB 
is very efficient, the values provide very small carbon (Xi2 C ~ 1.1 ■ 10~ 4 for 4-5M models 
down to Xi2 C ~ 1.6 • 10~ 5 at 6.3M ), and low oxygen (from X 016 ~ 6.3 • 10~ 4 at 4M to 
X016 — 10~ 4 at 6.3M ); nitrogen is very high, especially in the 4M , where the effect of the 
3DU is more visible (Xi4 N ~ 2.6 ■ 10~ 3 ). On the contrary, for the super-AGBs, we consider 
the values at the 2DU. The carbon abundance increases from Xi2 C ~ 8.9 • 10~ 5 for the 
7.5M model to Xi2 C ~ 8.1- 10~ 4 at 9M . The nitrogen is about constant (Xu N ~ 2.4-10" 4 ), 
whereas the oxygen abundance is constant (Xqiq ~ 3.8- 10~ 4 ) up to 8.5M , and increases up 
to Xqiq ~ 4.4 • 10~ 4 for the 9M model: this is due to the deeper sinking of the 2DU when 
increasing the mass. In the 9.5M model, the dredge-out process rises Xi2 C to ~ 1.8-10" 3 and 
X016 to ~ 7.7 • 10" 4 . The following evolution will reduce oxygen and carbon and increase 
nitrogen, if HBB is efficient in the se stars as it is in the mas sive AGBs (this depends on 
the convective model adopted; see IVentura fc D'Antonal 120051). As for t he 3DU, there are 
different and incomplete literature results (see Izzard fc Poelarends 20061 . for summary). In 
any case, it can only act to increase the values of total CNO abundances given in Figure [TJ 

How many stars would be affected if all the super-AGBs contribute to the second 
stellar generation? If all the super-AGBs contribute to the very helium rich population, the 
fraction of high-CNO stars, assuming a Salpeter's IMF, would be ~10%! This is a small 
fraction of the cluster stars, if we consider that the super-He rich stars are 15-20% of the 
total population, but spectroscopic observations of a large sample of stars could falsify this 
hypothesis. 



- 7- 



4. Do the most massive super-AGBs explode? 



Should the constancy of CNO be observationally confirmed also for the most helium 
rich population showing up in the quoted 4 clusters, we will be able to conclude that the 
most massive super-AGBs did not take part in the process of forming the second stellar 
generation in GCs. This may mean: 1) that there were reasons, independent from the 
super-AGB evolution, by which the self-enrichment process was effective either only at 
earlier ages (due to the envelopes of massive stars) or at a later age (during the normal, 
massive AGB evolution) 2) or that the self-enrichment was forbidden by the injection in the 
cluster gas of the energy due to the ecSN explosion of the stars which undergo the dredge 
out. Althou gh the explosion ene rgy of such event is significantly lower than inferred for core 
collapse SN (IDessart et al.l 120061 ). it is probably more than sufficient to expel the SN ejecta 
from most clusters. The conclusion then would be that at least a fraction of super-AGBs 
must evolve i nto ecSN, reinforcing the idea that ecSN channel for super -AGB stars could 



occur (see e.g. 



Siess fc Pumoll200a IPumo fc Siessl 12009 : IPoelarends et al.l 120070 . 



Neutron stars in Globular Clusters 



The final fate as ecSN of — at least a fraction — of super-AGBs is also required by the 
recent suggestion that NS in GCs should be mainly formed by the ecSN channel. 

The high number of neutron stars (NS) in GCs, whose presence is revealed mainly by 
the very high number of millisecond pulsars, remai ns a puzzle. The m assive GC 47Tuc was 
predicted to contain more than 1 000 NS (see e . g. IPfahl et al.l 120021 ) . Now this figure can 
perhaps be reduced to ~ 300-600 jHeinke et allboosh . The velocity distribution of young 
pulsars in the Galaxy shows that, at formatio n, the NS re ce ives a "nata l kick" , most likely 
due to the asymmetry in the supernova ejecta iFryerl (120041 ) . IPfahl et al.l (120021 ) have shown 
that only u p to 8% of NS can be reta ined in the clusters. The kick v elocity distribution which 
they used ([Hansen &: Phinneyill997l ) has been recently updated by iHobbs et al.l (120051 ). who 
find a higher mean vel ocity and velocity d ispersion. Using these new data, simulations of 
the retention factor by llvanova et al.l (120071 ) show that almost no NS can be retained for an 
escape velocity of 40Km-s _1 . A natural way to solve this problem, is to hypothesize that 
there i s a stellar populatio n which evolves into ecSN, that probably have much smaller natal 
kicks (llvanova et al.l 120071 ). possibly proportional to their reduced explosion energy. This is 
also suggested by the e xistence of a class of massive X-ray binaries with small eccentricities. 
van den Heuvell (120071 ) also supports the evidence of two channels for NS formation — with 
either small or large natal kicks — from an exam of the double neutron star binaries. 
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Fig. 1. — The triangles plus dashed line show the helium content in the ejecta of AGB 
stars as function of the C— O core mass, in the computations by P. Ventura & F. D'Antona 
(in preparation). The initial masses are, from left to right, 4, 4.5, 5, 5.5, 6, 6.3M . The 
triangles plus contin uous line sho w the helium abundance in the envelope after the 2DU in 
the super-AGBs by iPumol (120061 ). Masses from left to right are 7.5, 8, 8.5, 9, 9.5M . The 
circles show the ratio of total CNO abundance (Fcno) with respect to the initial value for 
the same models. 



